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Abstract. We provide a direct proof of a quadratic estimate that plays a central 
- - - ^ role in the determination of domains of square roots of elliptic operators and, as 

■ shown more recently, in some boundary value problems with boundary data. 

\ We develop the application to the Kato conjecture and to a Neumann problem. 

ps| ■ This quadratic estimate enjoys some equivalent forms in various settings. This 

gives new results in the functional calculus of Dirac type operators on forms. 
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1. Introduction 

The goal of this paper is first to present a self-contained and simple proof of the 
following quadratic estimate, and second, to convince the reader that this is a central 
estimate in this 

Theorem 1.1. Let n, m he positive integers, H. = L^(R", C™) and D, B he operators 
CN ■ on Ti satisfying the requirements (H). Then one has the quadratic estimate 

^ ■ r°° df 

(1) / \\t''BD{I + t^''BDBDy\\\^ - <\\uf, for all u en. 

^ One uses ( , ) and || || for the hermitian product and norm on Ti.. The hypotheses 

(H) consist of the following set of requirements. 

OO ! (HI) The operator D : D{D) — > is a homogeneous kth order differential oper- 

ator with constant coefficients. 
(H2) D is self-adjoint. 
^ ■ (H3) D is strictly accretive on its range, i.e. 

.c? ; \\V''u\\ <\\Du\\, for all u G D{D) n R{D). 

(H4) B is a. bounded operator on H. 

(H5) B is strictly accretive on R(-D): there is a constant 6 > such that 

Re{BDu,Du) >6\\Duf, for all ueD{D). 
(H6) (Off-diagonal decay) For every integer there exists Cat > such that 



O 



\u\ 



(2) Wt^BDil + t"'BDBD)~^u\\L2(^E) < C^(dist {E, F)/ty 

for all t > 0, whenever E,F C R" are closed sets, u eH satisfies supp u <Z F. 
We have set (x) := 1 -|- and dist {E, F) := inf {|x — ?/| : x e E,y E F}. 

In this paper, if A is a (densely defined) unbounded linear operator on Ti then 
□ (A), N(A), R(y4) denote respectively, its domain, null space and range. In (KB), 
V^u = {d°'Uj)^a\=k,i<j<m cousists of all the partial derivatives of u of order k. The 
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assumptions (H2,4,5) imply that BD has spectrum contained in a double sector of 
the complex plane centered around R and give boundedness of the operator in (H6) 
(See Proposition inH]). The constant in ([T]) depends on the implicit constants in (H). 

We mention right away that our interest is in operators B of multiplication by 
B{x), identified as a matrix having coefficients in L°°(R",C) , in which case (H5) 
is a form of Garding inequality. When D is first order, i.e. k = 1, and B is such 
a multiplication operator, then the off-diagonal decay (H6) holds true. Moreover, 
when k > 1, then (H6) is still satisfied in the case of most interest to us. (See Section 
O) However, we wanted to enlighten the observation that only (H6) is needed (in 
our arguments). We also stress that D is not assumed to be one-to-one. 

This theorem is proved in for first order D, i.e. /c = 1, as a corollary of another 
quadratic estimate. Our direct proof is shorter and simpler from the algebraic point 
of view, and also from the analysis point of view even though the same deep ideas 
are involved (Carleson measures, T{b) argument). Furthermore, our proof allows a 
simultaneous treatment of higher order D, i.e. k >2, which is new. 

The interest of proving a quadratic estimate is mainly in the following proposition 
corollary of results developed in [22] . 

Proposition 1.2. Assume that B,D satisfy (H2,4,5) on a Hilbert space Ti, that 
BD satisfies the quadratic estimate ([H), and that B*D satisfies the same quadratic 
estimate with B* in place of B. Then the operator sgn{BD) is bounded on Ti. and 
invertible on R{BD). 

The operator sgn^BD) is zero on n{BD) and satisfies {BBBDy/"^ = sgn{BD)BD 
on D{D). More is true, in particular BD has a bounded holomorphic functional 
calculus on Ti. We remark that the specific nature of Ti, B and D is not used in this 
proposition, which follows from operator theoretic considerations, once quadratic 
estimates for the operators BD and B*D have been proved. 

When k = 1, we obtain the following corollary to this result, once we have proved 
Proposition 15. 1[ Note that if B satisfies (H4,5), then so does B*. 

Corollary 1.3. Assume that B,D satisfy (Hl-5) on L^(R", C"), that k = 1, and 
that B is multiplication by a function B G L°°(R", £(C™)). Then the operator 
sgn{BD) is bounded on Ti. and invertible on R{BD). 

When A; > 1, we do not know if B being a multiplication operator is enough, in 
addition to (Hl-5), to conclude for the boundedness of sgn{BD). It is the case when 
B^D are as in Section [2TT] and the boundedness of sgn{BD) appears new. 

Known consequences of the boundedness of operators sgn(i?D) are short proofs of 
the Kato conjecture for elliptic systems in divergence form [HIS] and the boundedness 
of the Cauchy integral on Lipschitz curves [13] (see Section [2]). In Section [8] we give 
a pedestrian account of one of the results obtained in [3] concerning boundary value 
problems for second order elliptic systems which, in particular, give new proofs 
of solvability for single equations with real symmetric coefficients established in 
[ini HSl [20] . In Section [ini we show that the quadratic estimate for BD has different 
equivalent formulations with operators built by functional analytic considerations, 
including the one studied in [8], and we present a new application related to EVPs 
for differential forms. 

The quadratic estimate ([I]) has some further interest. It is easily seen to be 
stable under perturbation with lower order terms. This implies a simple proof of 
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the Kato conjecture for inhomogeneous elliptic operators (or systems) in divergence 
form, where previously it required an interpolation procedure from pure operators 
or a longer argument [Sj (See also [^). The extension to inhomogeneous situations 
is motivated also by potential applicability to time-harmonic Maxwell's equations. 
See the introduction of [3]. 
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2. Kato and Cauchy 

We present two typical applications of the boundedness of sgn{BD) already in 
the literature (at least when k = 1). We refer to [S] and the references therein for a 
number of further applications. 



2.1. Kato. The application to the square root of elliptic systems L = (V^YAV^ 
is as follows: A is multiplication by a bounded matrix A{x), and one assumes the 
Carding inequality 

Re{AV''u, V'^m) > 5|| V^wf , for all u G H^R", C^). 

Here u is C^-valued. Thus, we set ^^(R", C"*) = ^^(R", C^) ©^^(R", C^^) where 
m = N + Np and p is the length of the array V'^, 



(3) 



D :-- 



(V^)* 
V'^ 



B :-- 



I 
A 



One easily checks (Hl-5). For (H6), see Section El If M = A(V^')(V'')*, then 



{BDf 



L 
M 









M 



Since ^/{BD)^ = sgn{BD)BD, we get for u E L^(R",C^) under appropriate do- 
main assumptions that. 



Lu\ 



BD 



\\A{V'u)\\ 



" dx 



and 



2.2. Cauchy. As for the Cauchy integral, assume n = m = 1, D = 
B is multiplication by b{x) = where a G L°°(R, C) with Rea > 5 > 0. Then 
sgn{BD) is similar to the Cauchy integral on the Lipschiz curve with parametrization 
z{x) defined by z'{x) = a{x). 
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3. Proof of the main theorem 

3.1. Functional calculus for BD. First we need some review on functional cal- 
culus. Because of (H2), D is closed and densely defined and there is an orthogonal 
splitting 

(4) rL = n{D)®R{D). 
Define closed double sectors in the complex plane by 

S^:= {z eC ■.\±axgz\<ujU {0}}, 
and define the angle of accretivity of B to be 

uj := sup I a.Tg{Bv, v)\ < 7i/2. 

Proposition 3.1. Under (H2,4,5), we have 

(i) The operator BD is uj-hisectorial, i.e. a{BD) C 5*^^ and there are resolvent 
bounds II (A/ - BD)-^\\ < 1/dist (A, 5^) when A ^ S^. 

(ii) The operator BD has range R{BD) = BR{D) and null space N{BD) = N{D) 
such that topologically (hut in general non-orthogonally) one has 

n = R{BD) © N{BD). 

(iii) The restriction of BD to R{BD) is a closed and injective operator with dense 
range in R{BD), with estimates on spectrum and resolvents as in (i). 

These properties of closed operators of the form BD have been known for some 
time in the case when D is one-one, see for example p. When D is not one-one, 
first prove (ii), using (jlj) and (H5), and then adapt the proof in pQ to prove (iii). 
Part (i) follows. Note that this proposition only uses the fact that D is self-adjoint 
and B bounded and strictly accretive on R(-D). 

We set Rf = {I + isBD)-^ for seK. Then 

gf = ^(^^t^ - R^) = t^BDil + t^'^BDBD)"' 

and also \{Rt^u + Rfk) = {l + t^^BDBD)-^ . It follows from the previous result that 
i?f , hence Qf and (1 + t'^'^BDBD)^^, are uniformly bounded operators on 7i. 

We now come to the proof of Theorem 11.11 and assume all the requirements in 
(H). 

3.2. Reduction to a Carleson measure. Observe that by item (ii) of Proposition 
13. H as Qf vanishes on U{BD) it is enough to prove the quadratic estimate ([1]) for 
u G R{BD), hence for u € R{BD). Setting = QfB, it amounts to showing 

WOtDvf f < \\Dvf for all v G D{D). 
Let Pt be a nice approximation of the identity, i.e. the convolution with a real 
valued function t~"'(f{x/t) with (f smooth and having Fourier transform identically 
1 near 0. Let Pt act on C^-valued function componentwise. 

Proposition 3.2. 

/"°° dt 

(5) / \\Q,(l-P,)Dvf-<\\Dvf, veD{D). 

Jo 
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Proof. For the purpose of this proof, by using the sphtting (jl]), one can even assume 
V G R(-D)- Since Pf and D commute and (J — Pt)v G D(i5), we have 

Gtil - Pt)Dv = {QtD){I - Pt)v = t\BDf{I + {t'BD)Y\l - Pt)v. 
Now {t''BDf{I + {t^BDfy^ =1 -(1+ {t^BDfy^ is uniformly bounded, hence 

\\Q,{I-P,)Dv\\<^^\\{I-PMV 
Standard Fourier arguments show that 

and we conclude the proof of ^ using (H3). □ 

Remark 3.3. There are different possible choices of Pf's. For example, following 
[8] one can take Pt = {I + t^^D'^)~^. The organisation of the reduction to a Carleson 
measure would be somewhat different. 

Next, we perform the principal part approximation. 

We use the following dyadic decomposition of R". Let A = U^-oo where 
A2. := {T{k + (0, 1]") : k G Z*^}. For a dyadic cube Q G A2., denote by 1{Q) = T 
its sidelength, by \Q\ = 2"--?' its volume. We set At = Aaj if 2-?"^ <t <2K Let the 
dyadic averaging operator St : H ^ H he given by 

Stu{x) := Uq ■.=j u{y) dy= u{y) dy 

for every x G R" and t > 0, where Q is the unique dyadic cube in A^ that contains 
X. We remark that S"^ = St. 

Definition 3.4. By the principal part of {Qt)t>o we mean the multiplication oper- 
ators 7f defined by 

7f(x)w := {Qtw){x) 

for every w G C". We view w on the right-hand side of the above equation as the 
constant function valued in C" defined on R*^ by w{x) := w. We identify 7t(x) with 
the (possibly unbounded) multiplication operator 7^ : /(x) t— > 7t(x)/(x). 

Lemma 3.5. The operator Qt extends to a bounded operator from L°° into L|,^. In 
particular we have well defined functions '-ft G L|,^(R"; £(C"*, C™)) with bounds 

{bt{y)?dy<l 
JQ 

for all Q G A(. Moreover, ||7t5'(|| < 1 uniformly for all t > 0. 

Proof. Fix a cube Q G A^ and / G L°°(R™,C'") with ||/||oo = 1. Then write 
f = fo + fi + f2 + ■ ■ ■ where /o = / on 2Q and elsewhere and if j > 1, fj = f on 
2^~^^Q \ 2^Q and elsewhere. Then apply 6* and use (H6) for each term Qtfj with 
N large enough and sum to obtain 

h{Qtf){y)\'dy<C. 
JQ 
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If we do this for the constant functions with values describing an orthonormal basis 
of C™ and sum, we obtain an upper bound for the desired average of 7^. Next, for 
a function f eH, 



htsjr = E 



QeAt 



Q 



itiy)[ 4 f 



f 



< 



We have the following principal part approximation of by --^fSt 



Lemma 3.6. We have 



□ 



(6) 



\\QtPtf-ltStf\?j<\\f\\\ fen. 



Combining this with Proposition \3.2l we obtain the principal part approximation 



(7) 



r°° dt 

/ \\QtDv--itStDv\\''-<\\Dv\\\ veD{D). 
Jo * 



Proof. Write 

QtPt - itSt = {QtPt - itStPt) + iitStiPt - St)) + i^tS^ - jtSt). 

Because 5*^^ = St, the last term vanishes. Next, as •jtSt is uniformly bounded as an 
operator on H, we have 



f HStiPt - St)fr j<j^ \\{Pt - St)fr J 



< II fii^ 



The last inequality is done componentwise and is classical (See, e.g. [H], p. 172). 
We pass to the first term. We remark that for t > fixed and x G R", then 



{GtPt - ltStPt)f{x) = Qt[g -j^g^ {x) 
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where g = Ptf and Q is the only dyadic cube in containing x. Define Co{Q) = 2Q 
and CjiQ) = \ 2^Q if j G N*. Then 




<t'\\Vgr. 

We successively used the Minkowski inequality on the second line, (H6) on the third 
one, Cauchy-Schwarz on the fourth, Poincare inequality on the fifth, the covering 
inequality J^q^Ai ^ 2-'" and £{Q) ~ t on the sixth and the choice N > n + 2 

in the last. 
Hence 

\\e,pj - 7, wir j<f pvp./f j<\\fr 

using the standard Littlewood-Paley inequality on each component of /. □ 

Before we state the conclusion of this reduction, there is an essential observation. 
Identifying constant functions with their values, observe that Dv takes values in the 
vector space 

V = {DL : L: R" ^ C™, L a polynomial of degree k} C 

and so does St{Dv). Therefore, one considers the restriction of 7t(x) to V. Hence- 
forth, we consider 7t(x) as an element of C{V, C™) and its norm |7t(x)| is measured 
in this space. 

Recall that R"^""^ C™), (x, t) ^ 7t(a^), is a dyadic Carleson function if 

there exists C < oo such that 

// \it{x)?'^-^<C'\Q\ 

JJR{Q) ^ 

for all dyadic cubes Q C R". Here R{Q) := Q x (0, 1{Q)] is the Carleson box over 
Q. We define the dyadic Carleson norm ||7t||c to be the smallest constant C. The 
form of Carleson's lemma that we need and applied componentwise is as follows (see 
[6], p. 168 and references therein). 
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Proposition 3.7. 

htStDvrj<ht\\l\\Dv\\\ veD{D). 



Therefore, we have obtained 



Proposition 3.8. If the restriction of'-ft{x) to V is a dyadic Carleson function then 
the conclusion of Theorem \l.l\ holds. 

Remark 3.9. At this point, it is nowadays understood that the Carleson measure 
estimate can be achieved by what is called a T{h) argument, which consists in finding 
suitable test functions adapted to the operator 0t = QtB. However, a dichotomy 
appears on remarking that we can prove ([7]) for functions of the form / = Dv G R(-D), 
but not for functions / G N(D). This comes from the use of (H3) in Proposition I3.2[ 
The simple situation is when D is one-one (or, equivalently, D has dense range by 
(jlj)): the test functions are simply the columns oi B~^. When D fails to be one-one, 
this choice does not work as we have to select test functions in the range of D. For 
the Kato problem, one has one-oneness of the D involved only in one dimension. It 
is for this reason that the Kato problem was more difficult in dimensions n > 2 than 
in one dimension. For a fair comparison, we provide the concluding argument of the 
proof of Theorem 11.11 in both cases. 

3.3. The T(b) argument when D is injective with dense range. Fix a dyadic 
cube Q, and let tiq be a smooth real valued cutoff such that ?7q|2q = 1, supp (?7q) C 
3Q and || V^t/qHoo < for j = 1, 2, . . . , A; with / = 1{Q). Denote by Bj^ the j'th 
column vector in the matrix -B~^, and estimate 



l7.(x)r^< // \lt{x)S,{r^QB-^)\'^-^ 

R(Q) ^ JJR{Q) ^ 

J Jr(q) 



t 



< E 11(1 + t^\BDrrt'BD{r^Qe,)r^ + WvqBtY 

<E / \\D{riQe,)rt''-'dt+\Q\<\Q\. 



For the first row, we use the fact that B~^ is strictly accretive on 7i, and hence is 
pointwise uniformly strictly accretive. Here St acts componentwise on the matrix. In 
the first term of row four we write BrjqBj^ = rjqBBj^ = tiqCj, where Cj is the j'th 
standard basis vector in C"*. To obtain the second term, we apply the principal part 
approximation ([71), using the assumption that the range of D is dense in TC (hence 
Dv there can be replaced by any function in 7i). In row five we use the uniform 
boundedness of the operators (1 + t^'^(i?_D)^)~^ and that D{r]Qej) is supported on 
3(5 and is bounded by 
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3.4. The T(b) argument in the general case. We now consider the general 
case where D is not an injective operator with dense range in 7Y, so that we need 
to construct test functions which belong to the range of D. Fix Q a dyadic cube 
and w E V with \w\ = 1. Let L be a polynomial of degree k such that w = DL and 
sup3Q |(9"L(x)| < < |a| < — 1 and define wq := D(r]QL), where r]Q is the 

cutoff above. It follows that 

Wq e R(-D), wq\2q = w, suppwQ C 3Q and ||wq||oo < C. 

Next we define the test function bg ^ for e G (0, 1) by 

6^,, := Dv-^,, v-^^ := (/ + z{d)'BD)-\vQL). 

Lemma 3.10. There exists C > such that for each w E V with \w\ = 1, each 
dyadic cube Q C R" and each e G (0, 1), 



(8) 



(9) 



L\'' <CielY''\Q\ 



Q 



w\^<C\Q\, 



(10) 



w 



\^lx)S,blix)\^'^-^<Ce-^^\Q\. 



R{Q) 



t 



Proof. Using (/ + isBD)-^ -I = -isBD{I + isBD)'^ and r/gL G D{BD), we have 

v^^, - vqL = -t{d)\l + z{d)'^BD)-\BwQ). 

The properties of wq and r]Q and the boundedness of (/ + isBD)~^B imply ([H]). 
Applying D we get, 

b^e -^Q = -i{dfD{I + i{dfBD)-\BwQ). 

The properties of wq and the boundedness of sD{I + isBD)~^B imply ([9]). 

Next, let ip: R*^ [0, 1] be a smooth function which is 1 on (1 — t)Q, on Q'^ 
with II V'^(/^||oo < C{tl)~^ with t G (0, 1) to be chosen. We can write 



bT^.-w 



Jq Jq 



w) 



I + 11. 



Using and the properties of if together with Cauchy-Schwarz inequality, we 
obtain 

|J/| < CVi. 

For /, we can write using the properties ip and integration by parts. 



JR" JR" Jq 

where Dip is some £(0™", C'")-valued function bounded by C|| V^<y5||oo and supported 
in Q \ (1 — so that we obtain 

|/| < CeVt'^-^/^^ 
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Hence, choosing t = e, we have shown fllUp . 

Eventually, to prove fill I) , we can use the principal part approximation in Lemma 
13.61 (backwards) because bg^^ = Dvq^^ and H^q J| < 1 and it suffices to establish 

(12) // \e,b-g,{x)\'^-^<Ce~''\Q\. 

J Jr{q) ^ 

Now, 

Qtbg^, = t''BD{I + t^^BDBDY^BD{I + i{elf B DY\'qQL) 

= t^BD{I + t^^BDBD)-\l + i{elfBDy^{BwQ). 

= it/ elf {I + t^''BDBDy\elfBD{I + i{elf BDY^Bwq) 

Since {I + t^^ B D B D)~^ and {elf BD{I + i{elf BD^^ are bounded uniformly with 
respect to t and e/, we have 

||e,6^,J < C(t/elf. 

Integrating over t G (0, /] we obtain ffT21) . □ 

We now perform a sectorial decomposition and then a stopping-time argument to 
estimate the dyadic Carleson norm on 7t(a:). Cover £(I^, C"^) by a finite number of 
sectors C^,^ = {k e C{V, C"*) ; \k - \K\-f\ < z/|k|}, with 7 G £(©, C"), I7I = 1, and 
V G (0, 1). The number v is to be chosen later. Fix such a sector. It is enough to 
estimate the Carleson norm of 

7t(x) = l7i(x)ec^,.7t(a;). 

Pick w eV,w* E such that {'jw, w*) = 1 and \w\ = \w*\ = 1. For any k G C^^u, 
we have 

Re [|k|(7w,w*) — {k,w,w*)~\ < 

thus 

(1 — 1^)1^1 < Re{Kw,w*). 
Fix a cube Q. Applying this to 7t(x) with {x,t) G R{Q), we obtain 
(1 - u)\^t{x)\ < Re{^t{x)w,w*) 

< Re(7<(x)5i6^,,(x),ti;*) + mx)\ Re{j{w - Stb^^^{x)),w*) 

+ Re((7t(x) - mx)\i){w - Stbl,{x)), w*) 

< \ltix)Stbl,{x)\ + mx)\Rei^iw - Stb^^^ix)), w*) 

+ umx)\\w - S,bl,{x)\. 

Thus one needs smallness on Re(7(w — StbQ^^{x)),w*) and a control on the size of 
\w — StbQ ,,{x))\ on a large portion of R{Q). 

Lemma 3.11. There exists eo G (0, 1) such that for all e G (0, eg), any dyadic cube 
Q contains disjoint dyadic subcubes Qi with 

(13) $^|Q.|<(i-e)|Q|, 

i 

(14) Re(7(ti; - Stbl^,{x)), w*) < lOCv^, (x, t) G /2(Q) \ UR{Qi), 

(15) |«;-5t6^,(x)| < (a;,t) Gi?(g)\Ui?(g,)- 
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Here, C is the constant appearing in Lemma \3.1(A 
Assuming this, then we obtain 

(l-v-lOC^e-v^^ |7t(x)| < |7t(x)^t6^,,(x)|, {x,t) e R{Q)\UR{Q,). 

Choosing e and then v small enough (depending only on C, hence on (H)), we have 
shown for all Q with the corresponding Qi 

(16) Mx)\<2\^,{x)SrQMl {x,t)eRiQ)\UR{Q,). 

We finish with a classical observation: fix 5 > and 

Iff I ~ / \ i2 dxdt 

= sup / / \lt{x)\ —— < oo 

IV I J J{x,t)sR{Q'),t>S ^ 

where the supremum is taken over all dyadic subcubes of Q. Then, if Q' is such a 
cube and Q'- are the subcubes of Q' given by Lemma [3.111 

dxdt 



IL 



\lt{x)\' 



{x,t)eR{Q'),t>5 



t 



(17) J J{x,t)eR(Q'),t>5 1^ j J J(x,t)£R(Qr),t>s r 

KACe-^'lQ'l+AQY.m 

i 

<ACe~^''\Q'\+AQ{l-e)\Q'\. 

Hence, dividing by \Q'\ and taking the supremum over Q' we obtain Aq < 4Ce~^'^~^, 
and in particular 

1 ff |^^(3,)|2^<4c<e-2fc-i. 



\Q\J J{x,t)£R{Q),t>5 ^ 

This is independent of 5 > 0, hence we obtain the desired estimate by letting 5 tend 
to 0. 

It remains to prove Lemma 13.111 
Proof. We fix a dyadic cube Q. We assume e small. Observe that 



Re(7 {w-jhl)l,w*)<Ce'" 



and 

k-6nJ^ < C. 



/ 

JQ 



We subdivide dyadically Q and consider for the subcubes Q' both conditions 
(18) Re(7 [w -j bl)j , w*) > 10Ce'/\ 



(19) / 

JQ' 



w-blf>Ce~\ 



If one or the other holds, we stop and put Q' in the sought collection of stopping 
cubes {Qi). If none of the conditions hold, we subdivide Q' and iterate the test on 
subcubes. 
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We note that (x, t) e R{Q)\\jR{Qi) exactly means that w — StbQ^^{x) = w—f^, 6q 
for a non-stopping cube Q' . Thus fll4p and ffTSl) hold immediately. 

It remains to show ( fT3l) . Declare Qj of type 1 if ( fTSj) holds and of type 2 if ( fT9l) 
holds. We let Sj = XI where the sum is retricted to cubes of type j. We might 
count twice cubes of both types but that is not a problem. For cubes of type 2, we 
have 

J Qi J Q 

For cubes of type 1, we have 

lOCe'/' Si < ^ Re(7 ( ^ ^ - , w*) 

\JQ / \JQ\UQi / 

Using ffTUl) and the Cauchy-Schwarz inequality for the last term, we obtain 

where X = T,i/\Q\ G [0, 1]. The positive root of the corresponding equation is on 
the order of 1 — 81Ce for e small enough. Hence, 

X <1-Ce 

for e small enough. Thus, the total contribution of cubes of both types does not 
exceed {1 — Ce + £)\Q\, which gives ffTHl) (assuming C > 2 which we may). □ 

4. Historical comments 

The (almost) self-contained proof of Theorem 11.11 follows very closely the strategy 
of [1] which, of course, builds upon the ideas of many authors, and it also incorporates 
ideas from the various extensions of this argument found later on. Locating the 
origin of this and that can be subtle for the reader, so we devote this section to 
historical comments, giving appropriate credit for the crucial steps based on our 
understanding. We do not mention less recent progress and refer to P, H] for this. 

A strategy to solve the Kato conjecture in all dimensions was introduced and 
developed in [6] under kernel bound assumptions. It first involves the reduction to 
Carleson measures in such a context - which was named "principal part approxima- 
tion" in [8] - as described in Section 13.21 exploiting earlier ideas of Coifman- Meyer 
[T^ further elaborated in works of Christ- Journe and Semmes [T2l [23]. The present 
formulation of the principal part approximation is closer to the one in [8]. We have 
chosen this formulation for the simplicity of its proof (assuming minimal knowledge 
of Littlewood-Paley-Stein theory). The strategy of [Sj required the existence of a set 
of appropriate test functions in order to prove the Carleson bounds via the "T(b) 
theorem for square roots". In [6], Chapter 3 this existence was made an assumption 
called there "the class (5) assumption" . The construction of such a set was achieved 
for the first time in [18] to solve the two dimensional Kato problem. Our choice is 
close to this one, rather than the one used later in [Hj. However, we need to exploit 
the observation made in [2] that one can reduce the action of 7t(x) to a subspace, 
while this is not necessary in [18] or in [1]. 
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The importance of the inequahty (|TU|) . or at least an integral version of it, was 
pointed out in [H], Chapter 3, and found its roots in The kind of stopping- 

time argument providing an inequality like (|TU|) leading to fll7l) is developed for the 
first time in this context in [18]. It is mentioned in [1] that, in retrospect, this 
stopping-time argument is akin to an argument of Christ [H] devised for proving 
a local T{b) theorem for singular integrals. The conical decomposition done in the 
space of constants (C") to estimate the Carleson measure associated to '~ft{x) was 
the main new ingredient of [1^. This provided a means to build a different set of 
test functions to solve the Kato conjecture in all dimensions under kernel bound 
assumptions. The removal of such kernel bounds was achieved in [1], thus proving 
the Kato conjecture for second order operators in full generality. The idea which we 
use of doing the conical decomposition, not in the space of constants, but within the 
linear space of matrices to which 7t(x) belongs, is an important observation, made 
in [5] for proving the Kato conjecture for higher order operators and systems. 

Note that our argument is developed on R", while the one in [1] was pushed in [9] 
to Lipschitz domains for mixed boundary value problems. The case of Dirichlet and 
Neumann boundary conditions had been previously done in [7] by a direct reduction 
to the R" case. It would be of interest to adapt Theorem 11.11 to domains and to 
obtain new proofs and generalisations of the results just mentioned. 



5. Validity of off-diagonal estimates 

Proposition 5.1. If k = 1 then (H6) holds for all B and D with (HI, 2, 4, 5) when 
also B denotes multiplication by a matrix-valued function B G L°°(R", £(C™)). In 
fact, one even has exponential decay. 

The proof is inspired by the one in [1]. 

Proof It is enough to consider Rf = {I + itBOy^ for t ^ 0, as Qf = ^{Rf-R^t). 
Let d = dist {E,F). We have already proved uniform bounds. So it is enough to 
prove ([2]) under the assumption that \t\ < arf for some constant a > to be chosen. 
Assume u eH. with supp u G F. 
Write 

E ■={xeW : dist (x, E) < |dist (x, F)} 
and let (p : R" — > [0, 1] be a Lipschitz function such that suppc/? C E, (f\E = 1 and 

||Vv9||oo<4/d. 

Let 7] = e"'^'^/* — 1 and observe that 77 = on F and r] = e"*^/* — 1 > ^e""^/* on E, 
thus 

^e-'^/ii?fn|U.(^) < Wr^RfuW < || [r/ 1, i?f ]n|| 

using that r]u = 0. Next, 

[r] I, i?f ] = itRf [BD, r] I]i?f = itRfB[D, t] l]Rf 

and [-0,771] is multiplication by a function me""^'^/^ where m is supported on E with 
norm not exceeding Cad\\S/ip\\oo/i < 4Ca/t. Thus, using the boundedness of 
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\\r^R'^u\\<t\\[DMRtA 



<ACa{\\r^Rfu\\ + \\Rfu\\). 

Hence, choosing a small enough (independent of gives ||?7i?fM|| < H-RfwH < 
and this proves the proposition. □ 

Proposition 5.2. Let k >2 and B, D he as in Section IK7[ Then (H6) holds with 
exponential decay. 



Q 



B 



Ml 




_U2_ 





Proof. Observe that 

t''AV''{I + t'^''L)-^ui 

The off-diagonal bounds (El) for (/ + t^'^L)-^ and t^V'^il + t^'^L)'^ have been 
known for some time: see [16] where it is done for the semi-group e~* ^ instead 
of the resolvent. However, there is an argument using the spirit of the proof of 
Proposition 15. II working directly with Qf instead of Rt- From there the off-diagonal 
bounds for (/ + 1'^'' L)~^t^ (V^)* A follow from a duality argument changing A* to A. 
We leave details to the reader. □ 

6. Some functional consequences of the quadratic estimate 

6.1. Proof of Proposition 11.21 We refer to [1] for details on functional calculus 
for the class of operators under consideration here. Let us just say that there is a 
way of defining sgn{BD) using the following formula 

(20) sgn(i?Z^)/ = c r{t^BD)\l + t''BDBD)-'f ^ = c H {Qff f ^ 

Jo >^ Jo t 

with c"^ = u^^~^{l + u^^)"^ du. This comes from the fact that the function 
z I— > c {t^ zY [1 + t^^z^)"^ Y is holomorphic on C \ iR where it coincides with 
sgn(z), defined to be 1 on the right half-plane and -1 on the left half-plane, in other 
words, the holomorphic extension of the sgn function on the real line to C \ iR. 

By item (ii). Proposition 13.11 it is enough to define and prove boundedness of 
sgn{BD) on N{BD) and R{BD) separately. For / G N{BD) then Qf / = for each 
t, thus sgn{BD) = on H{BD). 

It is easy to see that the integral ( l20l) converges in norm in 7i for / G D{BD) fl 
R{BD), because then 

\\{t''BDf{l + t^''BDBD)-^f\\ < min(^^r'=). 



Since D{B D)r\R{B D) is dense in R{BD), this defines sgn{BD) on the latter provided 
one shows \\sgn{BD)f\\ < c\\f\\ for / G D{BD) H R{BD). 
Let / G D{BD),g G D{DB*). Then 

\{{QffL9)\ = \{Qf{Qff),{Qfr9)\ < WQffWWiQfYgW 

and applying the Cauchy-Schwarz inequality 
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The first factor is directly controlled by c||/|| by assumption. For the second factor, 
write 

{QfYg = (/ + t^^DB*DB*YH''DB*g = t^DB*{I + t^^DB*DB*y^g. 

We shall show in a moment that quadratic estimates for operators DB* are a con- 
sequence of the assumed quadratic estimates for B*D. We conclude that sgn{BD) 

is bounded as desired. 

We remark that sgn{B D)sgn{B D) = I on R{BD) from the properties of functional 
calculus. This gives the invertibility of sgn{BD) on R{BD), and the proposition is 
proved. 

6.2. Operators of type DB. 

Proposition 6.1. Under (H2,4,5), we have 

(i) The operator DB is uj-hisectorial, i.e. a{DB) C and there are resolvent 
bounds - DB)~^\\ < 1/dist (A, 5*^) when X ^ S^. 

(ii) The operator DB has range R{DB) = R{D) and null space N[DB) such that 
topologically (but in general non- orthogonally) one has 

n = R{DB) © N{DB). 

(iii) The restriction of DB to R[DB) is a closed and injective operator with dense 
range in R{DB), with estimates on spectrum and resolvents as in (i). 

(iv) // BD satisfies the quadratic estimate (tTp, for example if (Hl-6) are all 
satisfied, then for all g eH, 

/poo \ 1/2 

(21) n Wt'DBil + t^'DBDBr'grj] <\\g\\. 

If B is strictly accretive on all Ti, then DB = B^^{BD)B, so DB and BD are 
similar operators. In this case, bisectoriality, resolvent bounds, quadratic estimates 
and boundedness of functional calculus carries over from immediately BD to DB. 

However we are only assuming that B is strictly accretive on R{D) as in (H5). 
Denote by B' the restricted operator B' = B : R{D) BR{D) = R{BD). By 
(H4,5) this is an isomorphism, and we have DB = {B')~^{BD)B' on R{DB). Thus 
we can transfer results about BD on R{BD) to results about DB on R{DB). 

To extend the action of DB to all of Ti, we use (ii). To check (ii), note that DB 
and B*D are adjoint operators so that R{DB) = N{B*D)^ and N{DB) = R{B*D)^ 
and the splitting follows from Proposition 13. l( ii) on taking orthogonal complements. 

7. Spectral decomposition and dependence on B 

For this section, we specialise to the case of first order differential operators D, 
i.e. k = 1, and assume that S is a multiplication operator. 

We continue to treat DB, though the following results are readily adapted to BD. 

Proposition 7.1. Suppose the hypotheses (Hl-5) are satisfied with k = 1, and that 
B is multiplication by a function B G L°°(R"', £(0™)). Then the following hold. 

(i) The operator sgn{D B) is bounded on Ti. 

(ii) The operator DB has a spectral decomposition R{D) = R{DB) = T-lr)B+ © 
Hdb- where TCdb± = {v E R{D) : sgn{DB)v = ±v}. 
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(iii) The projections of R{D) onto 7ir>B± ore = 1/2(7 ± sgn{DB)) . 

(iv) The restriction of DB to 'Hr)B+ is uj-sectorial and thus generates a hounded 
analytic semigroup e"*-^^ in Hdb+ which satisfies lim^^oo and 
limt^o — > V for all v G T-Cdb+- 

Indeed V := e~*^^v is the unique function in C^(R+,7i) which satisfies 
§lV{t) + DBV{t) = for t > 0, as well as the limiting conditions just mentioned. 

Part (i) is proved in the same way as Proposition 11.31 once we have the quadratic 
estimates ([T]) and ( 12T1) . The other parts follow as a consequence. 

We remark that, by the preceding two propositions, there is a spectral decompo- 
sition 

n = nDB+ © Hob- © N(D5) 

corresponding to the right and left sectors of S^j and {0}. 

It is an important consequence of proving bounds for a general class of matrices B, 
that all of the corresponding bounded operators such as sgn{DB) depend analytically 
on B. To prove this, start by showing analytic dependence for the resolvents, and 
then prove analyticity for more general operators by using the fact that uniform 
limits of analytic functions are analytic. See [HI E] for further details. We need this 
fact for the projections defined above. 

Proposition 7.2. Suppose that z ^ B^ is an analytic function from z in an open 
set Q G C to functions Bz G L°°(R", £(C'")), and that the operators D and mul- 
tiplication by Bz in Ti satisfy (Hl-5) with k = 1 under uniform bounds. Then the 
projections E^" : R{D) Hdb^± C R{D) depend analytically on B. 



8. Elliptic systems with square integrable boundary data 

Let us illustrate the power of our estimates by showing that elliptic systems on 
R^+" = {{t,x);t> 0} of the form 

(22) diYt,,A{x)Vt,,F{t,x) = 

are well-posed under Neumann data on R", when the coefficient matrix A is 
self-adjoint, strictly accretive, and has coefficients which are bounded measurable 
functions of x G R". The functions F map R^^" to C^. One can also handle 
Dirichlet and Dirichlet data, thus generalising results concerning the case N = 1 
and real symmetric coefficients [IHl [151 ED]- See [3] for a more extensive theory and 
for the historical background. 



8.1. Results. On writing 
U{t,x) 



' Uo{t,x) ' 






U(t,x) 




_ V^F{t,x) _ 



Vt,xF{t,x) 



where U takes values in © C"^, the second order equation can be rewritten as 
a first order system 



(23) 
(24) 



diVt,xA{x)U{t,x) = and 

CUT\t,xU(t,x) = 
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for {t,x) e Rl+\ We assume throughout this section that the coefficient matrix 
A G L°°(R", C{C^ © C"^)) is strictly accretive in the sense that 





/ 




/ 









> K 



for some k > and all / e L2(R", C^), g G H\R'', C^). For the moment we do 
not assume self-adjointness of A. 

The Neumann problem for ( l22l) is well-posed in the L°°{L'^) sense means that, 
given w e L2(R'^, C^), there exists a unique function U G C^(R+, L2(R", © 
C"^)) which satisfies ^ and dMD on R^+", as well as limt^^U{t,.) = and 
limf^o U{t, .) = ti in the sense, where 



n 



(25) {Au)o 

Our aim is to prove the following result. 



j=0 



on R" 



Proposition 8.1. Assume, in addition to the above conditions, that A is self- 
adjoint. Then the Neumann problem for ^2M) is well-posed in the L°°{L?') sense. 

In fact the solutions satisfy quadratic estimates and have non-tangential maximal 
function estimates. See [3] for this, and for a treatment of more general conditions 
on A and other boundary conditions. 



8.2. A related equation. In what follows, we write 
(26) A -- 



a b 
c d 



where a, b, c, d are L°° functions taking values a{x) G £(C^), b{x) G £(0*^^, C^), c{x) G 
/:(C^,C"^),rf(x) G /:(C"^) for a.a. x G R'^. Also write 



(27) 



U 



u 



V 



V 



a b 
1 



U 



where f/o, Vq : R^^"" and U, V : R^+" ^ C"^ and set 



(2^ 



where V = Vt and div 



D 



div 
-V 



-V*. This operator D satisfies (HI, 2, 3). 
Note that R(D) = L2(R",C^ 



coercivity on R(-D), i.e. 



R(V). So the coercivity condition on A is really 
Re{ADu,Du) > ^pwlP 

for all u G □(£)). 

We note that this inequality implies the pointwise accretivity Re(a(x)^, ^) > 
for all ^ G so we may define 



'10" 




a b 


c d 




1 



In the next subsection we show that A is also bounded and strictly accretive on 
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Proposition 8.2. A function U G Ci(R+, /.^(R", © C"^)) satisfies ^ and 
g on R^+" if and only if the function V G C\li+, /.^(R", © C"^)) satisfies 



(29) 
(30) 

on R^+". 

Proof. Equation ( l29i) . namely 
d 



d 

-V + DAV 
ot 







and 



curLV(t, .) = 







V 


+ 





-V 



div 




dt 

is equivalent to 

d_ 

di 

whicli is in turn equivalent to the pair of equations 





"10" 




a b 


-1 


" Vo " 




c d 




1 




V 



a b 




" Uo " 


+ 


" div " 




"10" 




' Uo " 


1 




U 


-V 




c d 




U 



[| div] 



a b 




' Uo " 


a d 




U 







The first of these is fl2^ . while the second, when taken together with 
The Neumann condition (l25l) becomes Vq = on R". 

8.3. The map A 

strictly accretive matrix A 



IS 



□ 



A. Let us look at the correspondence between a bounded 

1 



before, the sub-matrix a is itse 
-1 



a b 
a d 



and A 



1 

c d 



a b 
1 



. As observed 



inverse 



a b 
1 



f bounded, strictly accretive on /.^(R", C^), and the 
exists and is bounded. For interest, we note that A = A. 

Proposition 8.3. Given a bounded matrix A as above which is strictly accretive on 
R{D), then A is also bounded, strictly accretive on R{D). 

Proof. This is easily verified, once we have the following identity: 



Re A 



a b 
1 



u 



a b 
1 



no 
u 



Re((uo, auo + bu) + {cuo + du, u)) 
Re((auo + bu, uo) + {cuo + du, u)) 





no 




Uo 




u 


5 


u 



□ 

Self-adjointness is not preserved under this transformation. Indeed, it can readily 

a P 



be checked that A is self-adjoint if and only if A has the form A 
a = a*, P = -7* and 6 = 6*. 



7 6 



with 
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8.4. Results for DA. We now apply the theory which we have developed concern- 
ing operators of the form DB to the present situation. So take Ti. = L^(R", © 
C"^), 5 = i and 

div 
-V 



D 



The pair [D, A) satisfies (Hl-5) and A is a multiplication operator. Therefore, 

H, 



by Proposition 17.11 the spectral projection : R(-D) 



'■DA+ 



is bounded. 



The functions v E 'Hjj^^ are precisely the boundary values of functions V on 11^"' 



which satisfy 
problem for 



and (1301) . Moreover V(t) = e ^^^v for t > 0. And the Neumann 
is the determination of and hence V , from its first component, 



^0 = where w G L^(R", C^) is given. We are writing v 

2CDn rtN 



V 



. So, on defining 



the operator : L (R", C ) by P^v = Vq, we have the following result. 

Proposition 8.4. The Neumann problem is well posed in the L°°{L'^) sense if 
and only if P^ ■ L^(R", C^) has a bounded inverse. 

This means that, given w G L^(R'^, C^), there exists a unique function v G 'H^^_^_ 
with Vq = w, and hence there exists a unique function V G C^(R+, L^(R'^, © 
C"^)) which satisfies |y + dAv = 0, \imt^^ V{t) = 0, lim^^o V{t) = v. 

Now the Neumann problem is not always well-posed [21], and P^ is not always 
an isomorphism. 

Our aim though is to prove Proposition 18. H or in other words, to show that when 
A is self-adjoint, then P^ is an isomorphism. 



8.5. Unperturbed operators. In order to tie our results in with the classical 

V 



theory, we start with A = A 

and Dv = ^/{D^)v, i.e. 



/. In this case v 



G 7iz)+ if and only if 



' div " 




^'0 




-V 




V 









V-Vdiv 



V 



where A = divV is the Laplacian on R". That is, v = — V(— A)^ ' i;o , or equiva- 
lently, Vq = (— A)~^/^divv and curlv = 0. 

Thus the map Pj is an isomorphism, and so, as we know, the Neumann problem 
for the Laplacian on the upper half space is well posed in the L°°{L'^) sense. 

8.6. Rellich inequality. We now consider bounded, strictly accretive, self-adjoint 

a (3 



matrices A. In this case A has the form A 



7 5 



with a = a*, (3 = —7* and 



(5 = 5*, or in other words, {A)*K = KA where K 
Since KD + DK = 0, we have, for v ETi 



'■DA+ 



1 

-1 
and V = e'^'^^v, that 



{KAv, v) 







or in other words, 



{K§:V,AV)+{KAV,i:V) 



{KDAV,AV)+{DKAV,AV) = 
(5v,v) = 0. 
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Therefore, as v G R(V) and 6 is strictly accretive on R(V), we have 

||vf < Re(5v, v) < \\vof + ||v||||t;o|| 
and hence the Relhch inequahty 

l|v|| < ll^^oll 

for all V G Ti.^^^. This in turn tells us that the mapping : v Vq satisfies the a 
priori estimates H-PifH ~ ||f || for all w G 

To prove surjectivity of P^, use the method of continuity (i.e. the invariance of 
semi-Fredholm index) for the family Br = tA + (J — r)/, < r < 1, and the fact 
that Pi is an isomorphism. For this we need to know that the spaces 'hiDBr+j 
other words the projections depend continuously on r. 

In fact Bt- depends analytically on r G C. Moreover, on some open subset C C 
which contains the closed real interval [0, 1], the operators B^ satisfy (Hl-5) with 
uniform bounds. Therefore, by Proposition 17.21 the projections E^'" are actually 
analytic in r. 

This completes the proof of Proposition 18.11 

Remark 8.5. The results on analytic dependence imply even more, namely that 
the Neumann problem is well-posed in the L°°{L^) sense for all coefficient matrices 
A which are sufficiently small perturbations of a self-adjoint matrix. See [3]. 

Our aim has been to show that the functional calculus provides a mechanism 
to turn a Rellich estimate (i.e., comparability in the norm of transverse and 
tangential derivatives, which always holds for self-adjoint, t- independent, divergence 
form elliptic operators) into solvability. Previously, such a mechanism had existed 
in the case of a single equation with real symmetric coefficients (e.g., as in p!9l [20]). 
or in the case of a block diagonal matrix (i.e., the setting of the Kato problem, in 
which case the square function estimates are equivalent to a Rellich estimate), or in 
the case of a system with constant coefficients. 



9. Inhomogeneous D when k = 1 

Consider operators D and B on Ti. = L^(R",C"') with the requirement (inhH) 
which consists of (inhHl): The operator D : D{D) ^ 7i is a (inhomogeneous) first 
order differential operator with constant coefficients, (H2), (inhH3): 

Ml + ||Vm|| < \\Dul for all u G D{D) n R{D), 

(H4,5) and 5 is a multiplication by a matrix-valued function. 

Then we claim that the conclusion of the main theorem is valid replacing (H) by 
(inhH). The proof consists in going back to the homogeneous case (It would be nice 
to have a direct proof as in Section [3]). 

Write D = Di + Dq where Di is a homogeneous first order differential operator 
and Dq is multiplication by a constant matrix. We observe that since D is self- 
adjoint, so are Do and Di as seen by computing the Fourier symbols. Also, using 
Fourier arguments again, one can check that (H3) holds for Di. So one can define 
(/ + itBDi)~^, the corresponding Qj, and obtain the conclusion of Theorem 11.11 for 
Ql- Observe that 

\\{I + itBD)~^ - {I + 'itBDi)-^\\ = + itBD)-\tBDo){I + ztBDi)-^\\ < \t\. 
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Recall that = ^ii^-t - Rf)^ hence 

\\tBD{I + t^BDED)-^ - tBDi{I + t^BDiBDiY^\\ < t. 

So /o lletwf f < llwf for all M e 7^ follows from ([I]) for 

It remains to control ||6tM||^ y and we know that it suffices to assume u = Dv 
with V G 0{D). We can also assume v G R(-D) from (jlj). But ||6(£)|| < hence 
J-j°° ||GtDw||^ y < ||f P and we conclude invoking (inhH3). 

For the application to the Kato problem for second order elliptic operators with 



divAV + div6-cV- 
A 



lower order terms L = 
erators with L°°(R", C) coefficients and 
we take 



-d, where A, b, c, d are multiplication op- 
^ is strictly accretive on L^(R", C""^^), 










div^ 


l" 




'l 








(31) 


D := 










, B:= 





A 


b 






1 













c 


d 



Remark 9.1. The inhomogenous version for higher order Z) is in [6], Section 0.7]. 
Further results on inhomogeneous problems are in [^. 

10. Related operators and further applications 

The aim of this section is to see that operators BD or DB are the building blocks 
of other operators obtained by functional analytic considerations in such a way that 
results for DB and BD apply immediately. We will finish with a concrete application 
to the functional calculus of Dirac type operators on forms. 

10.1. Operators of type F + B~^r*B. This class of perturbed Dirac operators 
was studied in [8], where bisectoriality and quadratic estimates where proved. As a 
corollary, our Theorem 11.11 was proved in [HI Theorem 3.1(iii)]. In this section we 
prove the converse to this, i.e. we deduce the results in [H] from Theorem ll.il 

As in [H], we consider the following slightly more general situation. Let F and F* 
be adjoint nilpotent operators in an arbitrary Hilbert space 7i, i.e. F^ = (F*)^ = 0, 
let Bi, B2 be bounded operators such that Bi is strictly accretive on R(F*) and B2 is 
strictly accretive on R(F). Furthermore, assume that F*i?2-BiF* = and Fi?ii?2F = 
0. This holds in particular if {Bi, B2) = {B^^, B) with B strictly accretive on all of 
Ti, . For more details on this hypothesis, we refer to (Hl-3) in [8]. In this case, the 
operator on H 

Hb-F + F^, 
induces a Hodge type splitting 

(32) 7^ = (N(F^) n N(F)) 

Now consider the operators 

-Q Y* 



B 



B^T*B2 



RiT%) © R(F). 



D 



F 



B :-- 








B2 



acting in © 7Y. Note that the hypothesis on F, F*, B2 above is equivalent with 
(H2,4,5) for D, B. Using the fact that the Hodge splitting (l32l) is topological, the 
map 

" fi " 



S:n^n®n:fo + fi + f2 



-B2/2 
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is well-defined and bounded, and one can show that the restriction of S to R(n5) = 
R(r^) © R(r) is an isomorphism onto R{BD). Since 



BDS{0 + /i + h) = 



' fl ' 






-B2/2 




.i?2(r/i)_ 



= s{o + ry, + r/i) = siUsio + a + A)), 



we have shown that the restrictions BD : R{BD) R{BD) and Ub '■ R(IIb) — ^ 
R(IIb) are similar operators. Consequently, we obtain the following. 

Proposition 10.1. Let Ub in Ti and BD in Ti (B'H he as above. Ub is a bisec- 
torial operator on H with resolvent bounds. Furthermore, if BD satisfies quadratic 
estimates, so does Ub. 

10.2. Operators of type BDi + D2B^^. In this section we aim to deduce quadratic 
estimates for operators well adapted to boundary value problems for differential 
forms. Similar to our discussion of the class of operators 11^ above, we consider the 
following slightly more general class of operators 

B^Di + D2B2 

acting in an arbitrary Hilbert space Ti. We assume that the two pairs of operators 
{Di,Bi) and {D2, B2) both satisfy (H2,4,5), and the following compatibility condi- 
tions. For the unperturbed operators Di and D2, we assume that R(-Di) C N{D2) 
and R(-D2) C N(Di). This means in particular that D2D1 = = D1D2 on appropri- 
ate domains. For the perturbed operators we assume that B2B1 : R(-Di) — > N(Z}2) 
so that D2B2B1D1 = on D{Di). This holds in particular if {Bi, B2) = {B,B-^) 
with B strictly accretive on all of Ti. 

Similar to the fact that an operator acts by swapping the two ranges in its 
Hodge splitting, an operator B1D1 + D2B2 acts diagonally in the associated splitting 

(33) n = (N(fiiDi) n N(D252)) © R(fiiDi) © R(d^ 

of Ti. To see this splitting, note that 

n = R(5iDi) © N(5iDi) = U{D2B2) © R(D^, 

and (1331) follows by intersecting these two splittings since R{BiDi) C N{D2B2) and 
KD2B2) C N(fiiDi). Since BiDi and ^252 act as 

BiDi = © BiDi © 0, ^2^2 = © © D2B2 

in the splitting ( !33|) . the operator BiDi + D2B2 is the direct sum of these two 
operators, namely 

B^D, + ^2^2 ■.n^H:fo + h + f2^Q + B.D^h + ^252/2. 
This shows the following. 



Proposition 10.2. Assume that Di, Bi and D2, B2 are as above. Then B1D1+D2B2 
is a bisectorial operator on Ti with resolvent bounds. Furthermore, if BiDi and D2B2 
satisfy quadratic estimates, then so does BiDi + D2B2. 
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10.3. An application to Dirac type equations. We end with an application of 
the quadratic estimates for the class of operators BiDi + D2B2 to boundary value 
problems for differential forms. The goal is to prove that the underlying operator 
Tb used in [2] to obtain perturbation results for BVPs for Dirac type equations, and 
in particular for BVPs for second order divergence form equations, has quadratic 
estimates for all complex strictly accretive coefficients B. In [2], quadratic estimates 
were proved only for special types of coefficients B, namely for small perturbations 
of real symmetric, constant and block form coefficients. 

The operators Tg are infinitesimal generators for the studied Dirac equations, and 
following [21 Definition 2.10] we have 



B 



-iM^\d + B-^d*B), 



which acts in L^(R"'; A), where A = AcR-^^" is the full complexified exterior algebra 
of R^+". Here dis a nilpotent differential operator, i.e. d^ = 0, and 5 is a strictly 
accretive multiplication operator so that the operator Tg, modulo the factor —iMs 
(being an invertible, non-accretive, multiplication operator), is of type Ub- More 
precisely, if /i/ = cq a / is exterior multiplication by the basis vector normal to R'^ 
and m := fi + fi* then d := imd where d is the exterior derivative, and Mb '■ = 
fi*fi — B~^fifi*B. In the work [2], this factor complicated the application of results 
for the class 11^ and only gave partial results. However, in connection with the later 
work [3], it was realized that the operators Tg actually are similar to operators of 

type BiDi + D2B2. This similarity uses the transform B ^ B := BB of strictly 
accretive matrices analogous to Proposition 18.31 on splitting the space L^(R";A) 
into normal and tangential forms and writing the operators B, B_, B as the matrices 



B 



Bj_j_ B^ 

Bn I Bn 



B 



B„ 





-Bill 



B 



B , , B , 



This is summarized in the following new result. 



Corollary 10.3. Let B G L°°(R"; C{A)) be any complex coefficient matrix function 
which is strictly accretive on L^(R", A), matrix function. Then 



B 



B-' 



(^D,B + B-'D2^B, 



where the differential operators are Di := fi*d — fid* and D2 := fi*d* — f^id. In 
particular Tb is an injective uo-hisectorial operator, uj being the angle of accretivity 
of B, has resolvent bounds and satisfies quadratic estimates in L^(R"; A). 

Proof. Since Di,D2 are first order differential operators, according to Theorem ll.il 
Proposition 110.21 and Proposition 15. 11 it suffices to verify the similarity. Multiplying 
the equation with BMb, it suffices to show that 

-i{Bd + d*B) = {BMb)(B~^DiB + B-^D2B). 
Identifying cq a /i + /s G ^^(R"; ARi+") with [/i, /s]* E ^^(R"; AR")^, the above 



definitions give —id 
d' 

d ^^^^^ 



d 
-d 
^0 d 
d* 



—id* 



d* 
-d* 



BM, 



B 



-B, 






5„i 



. The similarity is now straightforward to verify. □ 
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We end with the remark that it is only the part Tb similar to DiB of the full op- 
erator Tb that is needed for the application to boundary value problems for /c-vector 
fields / /c- forms. This application is described in [31 Section 6]. The complementary 
part Tb similar to B~^D2 will mix /c-vector fields of different order k, but has the 
advantage of making the full operator Tb injective. 
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